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Simulation  and  Interpretation  of  Polarization 
Diversity  Radar  Spectral  Functions 


I.  INTRODUCTION 


The  research  described  in  the  following  sections  is  a  continuation  of  the  re- 

1  2  3 

search  reported  by  Metcalf  and  Echard  and  Metcalf.  '  The  earlier  work  docu¬ 
mented  various  theoretical  and  practical  aspects  of  the  interpretation  of  signals 
receivable  by  coherent  polarization-diversity  radars,  including  the  formulation  of 
the  received  signals  and  the  spectral  functions  derivable  from  them  and  the  effects 
of  propagation  phenomena,  non-Rayleigh  scattering,  variations  in  apparent  or  actual 
shape  and  orientation  of  hydrometeors,  and  air  velocity  variance  on  the  received 
signals  and  spectral  functions. 

To  elucidate  these  effects  we  developed  numerical  simulations  of  the  spectral 

4 

functions,  following  the  approach  of  Warner  and  Rogers,  who  computed  power 


(Received  for  publication  26  April  1983) 

1.  Metcalf,  and  Echard,  J.D.  ( 197  8)  Coherent  polarization -diversity  radar 

techniques  in  meteorology,  J .  Atmos.  Sci.  ^5:2010-2019. 

2.  Metcalf,  J.T.  (1981)  Propagation  effects  on  a  coherent  polarization-diversity 

radar.  Radio  Sci.  16:1373-1383. 

-  vW 

3.  Metcalf,  J.  1.  ( 1 983)  Interpretation  of  simulated  polarization  diversity  radar 

spectral  functions,  Radio  Sci.  IThin  press. 

4.  Warner,  C. ,  and  Rogers,  R.  R.  (1977)  Polarization -diversity  Radar;  Two 

Theoretical  Studies,  Sci.  Reprt.  MW-()0,  Stormy  Weather  Group, 

McGill  I'mversitv,  Montreal,  Quebec,  Canada. 


7 


spectra  corresponding  to  those  which  would  be  observed  in  the  in  the  transmission 
channel  (S^)  and  orthogonal  channel  (S2 )  of  a  10-cm  radar  transmitting  circular 
polarization.  The  present  model  yields  the  two  power  spectra  and  the  cross - 
spectrum  Sjg  of  the  two  received  signals.  The  latter  function  is  of  significance 
when  there  is  either  non-Rayleigh  scattering  or  a  non -negligible  propagation  effect, 
both  of  which  are  typically  encountered  at  shorter  radar  wavelengths.  The  model 
development  was  originally  conceived  as  a  counterpart  to  empirical  studies  based 
on  data  from  the  8.  6 -mm  coherent  polarization -diversity  radar  operated  by  the 
Wave  Propagation  Laboratory  of  the  National  Oceanic  and  Atmospheric  Administra¬ 
tion.  5  The  simulation  work  also  permitted  us  to  develop  and  refine  procedures  for 
manipulating  and  displaying  ;he  spectral  functions  which  could  be  used  with  spectral 
functions  derived  from  radar  data.  This  report  documents  the  formulation  of  the 
model  and  presents  a  variety  of  computational  results. 

g 

A  similar  set  of  model  computations  was  presented  bv  Metcalf  and  Matthews. 
However,  there  were  two  errors  in  the  model  formulation  used  in  that  report  which 
necessitated  the  publication  of  the  present  results.  One  error  was  in  the  computed 
magnitude  of  the  propagation  term,  affecting  computations  for  non -zero  propagation 
distance.  The  other  error  was  a  misinterpretation  of  the  scattering  amplitude  coef 
ficients  which  resulted  in  an  incorrect  orientation  of  the  scattering  amplitude  ratio 
in  the  complex  plane.  The  computed  cross  spectrum  was  thus  in  error  in  all  cases 
and  the  power  spectrum  in  the  transmission  channel  (S j )  was  in  error  in  cases  of 
non -zero  propagation  distance.  Most  of  the  qualitative  conclusions  of  the  previous 
report  remain  valid,  but  quantitative  details  differ,  particularly  in  cases  involving 
differential  propagation  effects.  In  addition  to  correcting  the  earlier  results,  the 
present  report  contains  results  of  computations  for  7.5  cm  wavelength,  correspond¬ 
ing  to  predominant  Rayleigh  scattering. 

The  results  shown  in  Section  3  illustrate  the  effects  of  air  velocity  v  ariance, 
noise,  and  differential  propagation  parameters  on  the  spectral  functions  at  a  par¬ 
ticular  elevation  angle  and  two  rainfall  rates.  These  results  emphasize  the  require¬ 
ment  foi  low  noise  in  the  received  signals,  so  that  the  spectral  power  ratio  Sj/S9 
and  the  cross -spectral  amplitude  ratio  S^/S^  can  provide  useful  information  on 
backscatter  characteristics  across  the  Doppler  velocity  domain.  The  presence  of 
noise  precludes  the  measurement  of  either  of  these  ratios  at  very  low  Doppler 
\elocities  where  propagation  effects,  if  present,  would  first  dominate  backscattering 


5.  Pasqualucci,  F.  ,  Bartram,  B.  W. ,  Kropfli,  R.  A. ,  and  Moninger,  W.  R,  (1983) 

A  millimeter-wavelength  dual-polarization  Doppler  radar  for  cloud  and  pre¬ 
cipitation  studies,  J,  Clim,  Appl.  Meteor.  22}758-765. 

6.  Metcalf,  J.  I. ,  and  Matthews,  J.  E.  ,  Ill  (1981)  Numerical  Simulation  of 

Coherent  Polarization-diversity  Radar  Spectral  Functions.  Scientific 
Report  No.  1,  Natl.  Sci.  Foundation  Grant  ATM-8018382,  Project  A2818, 
Engr.  Expt.  Stn.  ,  Ga.  Inst,  of  Tech. ,  Atlanta,  Georgia. 
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effects.  Hence  the  separation  of  propagation  and  backscattering  effects,  as  postu- 
lated  by  Metcalf,  will  prove  to  be  a  more  complicated  process  than  initially 
envisioned.  With  increasing  air  velocity  variance,  the  numerical  characteristics 
of  the  spectral  functions  at  all  velocities  tend  to  approach  the  characteristics  of 
the  corresponding  integrated  quantities.  Since  the  first  maximum  and  minimum 
of  the  Mie  scatter  cross-section  occur  at  relatively  large  drop  sizes  at  8.  6  mm 
radar  wavelength,  the  spectral  features  associated  with  these  extrema  occur  at 
Doppler  velocities  well  above  those  of  the  spectral  peaks;  with  increasing  air 
velocity  variance,  these  features  appear  at  relatively  higher  Doppler  velocities. 
Therefore,  in  the  interpretation  of  spectra  computed  from  experimental  data,  it 
will  be  necessary  to  make  an  adjustment  for  air  velocity  variance  before  particular 
features  of  the  spectral  functions  can  be  identified  with  particular  Doppler  fall 
speed  components. 


2.  NUMERICAL  CONCEPTS 


2. 1  Spectral  Function  Formulations 

The  power  spectra  derived  from  signals  in  the  transmission  channel  (Sj)  and 
the  orthogonal  channel  (S2 )  and  the  cross -spectrum  (Sj2>  of  the  two  received  signals 
can  be  expressed  as 

OO 

S.(v)  f  a  (r.  0)  G(V(r)  sin  0,  E  )  N(r) 

1  0 


x  |  v(r,  0)  eJl«<r.0)±2a<r,0)]  +2pej<\±  2t)|2  dr 
00 

S9(v)  /  a(r,  0)  G(V(r)  sin  0,  E  )  N(r)  dr 
0 

or 

S,9(v)  =  f  ff(r,  0)  G(V(r)  sin  O,  S  )  N(r) 

1  0 


(1) 

(2) 


(„<r.®>  ei[«<r.*>±2a<r.*>]  +  2pe->(x±2T>)  dr 


(3) 


where  v  is  the  Doppler  velocity  defined  as  positive  toward  the  radar,  a  is  the  back- 
scatter  cross-section,  r  is  the  equivalent  spherical  radius  of  a  drop,  0  is  the  eleva¬ 
tion  angle,  N(r)  is  the  drop  size  distribution  function,  and  a  is  the  orientation  angle 
or  canting  angle  of  each  scatterer.  The  velocity  weighting  function  is  defined  by 
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G(V(r)  sinp,  2)  = 


1 

-  exp 

Zy/TH 


( 


(v  -  V(r)  sin  0r 
2  22 


(4) 


where  V(r)  is  the  raindrop  fall  speed  in  calm  air  and  2"  is  the  Doppler  velocity 

variance  due  to  turbulence  or  wind  shear.  One  could  include  a  non -zero  mean  wind 

bv  substituting  [V(r)  sin  o  +  V  ]  for  [V(r)  sin  oj  in  Eq.  (4),  where  V  is  the  radial 

a  a 

component  of  the  air  velocity.  The  complex  backseatter  amplitude  ratio  is  defined 
for  a  spheroidal  hvdrometeor  bv 


e(r,  o) 


i 6(r, 0) 


S  (r.0)  -S  (r,  o) 

XX  r _ vv  _ 

S  (r.  0 )  +  S  (r,  oT 
xx  vy 


(5) 


where  S  ...  and  S  arc  the  b lekscattor  amplitude  coefficients  for  signals  of  polariza¬ 
tions  parallel  and  perpendicular,  respectively,  to  the  svmmetrv  axis  of  the  spheroid. 

For  Ravlcieh  scattering  from  dilate  spheroids,  S  >S  and  the  amplitude  ratio 

w  xx 

lies  near  the  negative  real  axis.  '!  r.is  con.at  ior  can  be  deseribi  1  «  ttiier  bv  v  <  0  and 
tj  -  0  or  bv  r  >  0  .,nd  5  -  ;Ve  use  tif  latter  1  onventi  >n.  The  propagation  term 
'  T*  is  base  i  on  a  m ridel  of  an  anisotropic  ( nreeinitatio'  filled!  medium  with 
svmmetrv  ixos  t  ot;:»--l  through  an  angle  r  from  the  vertical  and  horizontal.  The 
other  propagation  parameters  are  related  to  the  total  one -wav  differential  attenua¬ 
tion  A  A  an.t  lift''  r-mtial  phase  shift  .\<J>in  the  medium  to  firs'  order  approximation 
bv 


in1''  -  tann(0.  05? 50  A  A  :  j  .160)  .i'll)  .  *6) 

I  he  upper  in  i  lo  ver  signs  in  E  js.  (1)  and  CD  and  elsewhere  correspond  to  the  trans¬ 
mission  of  gin  am  left  i  treular  polarizat ion,  respectively. 

I  he  sp«  :  rum  defined  in  Eq.  <2>  vields  the  relation 

nr 

»J  /  S,,(v  t  i Iv  (7) 

-or  “ 

where  if  is  the  refleetivilv  i.  units  of  inverse  length,  that  is,  cross -section  per 
eiit  volume.  The  backseatter  er  iss-seetion  is  given  bv 

a  -  -  ,S  (r,  o)  ♦  S  (r.  o) 1  “  .  (ft ) 

xx  vy 
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For  our  computations,  we  used  values  of  S  and  S  computed  by  Oguchi  and 
7  xx  yy  h 

Hosova.  These  computations  were  based  on  the  linear  size -shape  relation  of 

8  9 

Spilhaus,  rather  than  the  more  realistic  Pruppacher  and  Fitter  size -shape  rela¬ 
tion,  and  our  results  therefore  cannot  be  strictly  representative  of  what  one  might 
observe  in  nature. 

4 

Warner  and  Rogers  computed  power  spectra  up  to  proportionality  for  Rayleigh 
scattering  from  raindrops  having  shapes  as  described  bv  Pruppacher  and  Pitter. 
Their  computations  included  air  velocity  standard  deviation,  radar  elevation  angle, 
rainfall  rate,  and  fraction  of  preferentially  oriented  scatterers  as  computational 
parameters.  The  elevation  angle  affects  the  relative  significance  of  fall  speed 
\  ariance  and  air  velocity  variance,  and  the  rainfall  rate  affects  the  drop  size  dis¬ 
tribution  through  the  assumed  Marshall  and  Palmer^  distribution  (unction.  We 
included  all  these  parameters  in  our  model  and  added  propagation  distance  as  a 
Parameter  to  simulate  non -negligible  differential  propagation  effects. 

The  functions  of  greatest  interest  arc  the  power  spectrum  (in  the  channel 
•  pposito  to  the  transmission  channel),  the  spectral  power  ratio  S./S  ,  the 
•tore. alized  eross-spei  truni  S  9/(S  .S9)  “,  and  the  function  S  which  is  of  par- 

■i  ttlar  interest  where  propagation  effects  are  present.  1'he  functions  S./S0  and 
Sj.,  S.,  arc  essentially  weighted  averages  of  the  squared  absolute  value  in  Eq.  (1) 
an  i  the  qu.r  :itv  in  brackets  in  Eq.  (3),  respectively,  and  can  be  expressed  as 


Sj  S.  r”  ■  -4(r  ■  4  If 


e  (  pe'|(x  1  2T)  . 1(5  •  2lY>  ) 


S,.,/S.,  -  r  e'1  {bl2‘y)  ’  2pcl<X  :2  Ti 


(9) 

(10) 


where  the  iverbars  indicate  integration  over  drop  size,  as  shown  explicitly  in  Eqs. 
( 1 ),  (2i,  and  (3).  Because  the  function  contains  the  backscatter  amplitude 

ratio  i  and  the  propagation  parameter  p,  which  is  also  an  amplitude  ratio,  we  sug¬ 
gest  that  this  function  he  called  the  "cross -spent ra!  amplitude  ratio.  "'  (The  com - 

larwble  parameter  for  non -coherent  radar,  W/W„  in  the  notation  of  McCormick 

11  ,  i  ,, 
and  llendrv,  ran  be  called  the  cross -covariance  amplitude  ratio.  )  The  nor¬ 
malize  i  cross -spectrum,  the  magnitude  of  which  is  the  coherency,  is  related  to 
tite  two  functions  above  bv 


o  /,*  c  .1/2  .  S12/S2 

b12/(t>lV  ‘  ~~7~~TT? 


(Sj/S,,)1 


(11) 


(Due  to  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  References,  page  51.) 
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For  modeling  of  the  spectral  functions,  it  is  useful  to  express  the  averaged  quan¬ 
tities  in  terms  of  parameters  that  allow  some  separation  of  variables  and  that 
characterize  both  the  extent  to  which  the  medium  is  preferentially  oriented  and  the 
orientation  angle.  Ultimately,  it  will  be  necessary  to  include  the  effect  of  shape 
variation  on  the  backscatter  power  and  amplitude  ratios.  The  following  subsections 
discuss  in  more  detail  the  formulation  of  these  ratios  in  the  model. 


2.2  Backscatter  Amplitude  Ratio 

The  backscatter  amplitude  ratio  defined  in  Eq.  (5)  is  shown  in  Figure  1  for  two 

wavelengths  and  two  elevation  angles,  based  on  the  Oguchi  and  Hosoya  amplitude 

coefficients.  Comparing  the  curves  for  the  two  elevation  angles  at  each  wavelength, 

2 

one  can  verify  that  the  magnitude  of  the  ratio  is  approximately  proportional  to  cos  0 
and  that  the  phase  angle  6  is  not  strongly  dependent  on  the  elevation  angle,  except 
for  drop  diameters  near  half  the  wavelength,  that  is,  between  4.0  and  4.5  mm 
diameter  for  8.6-mm  wavelength. 

To  evaluate  the  term  i‘  e-1^  1  ^  01  *  for  a  collection  of  raindrops  we  use  the  ap¬ 
proach  of  McCormick  and  Hendrv.  11  We  assume  that  drops  of  a  given  radius  have 
a  distribution  of  canting  angles  and  integrate  the  quantity  t'.fr,  0)  exp [ j(6.(r,  0)  ±  2  cy. (r,  ©))], 
corresponding  to  an  individual  drop,  over  all  orientations  to  obtain  the  result 


i' 


e,|(S  :  2  <V> 


j(6.  (r.  pi  i  2  tv(r,  ©)) 
p  (r.  o)  v.  (r,  0)  e  1 

<»  l 


(12) 


where 


-12 


p  (r.o)  -  f  cos  ( 2  (nr-  <v) I  T(a-d)  ci(a- or) 

,y  r. !  2 


(13) 


<v  is  the  average  orientation  angle,  and  T(<y-<y)  is  the  symmetrical  distribution  of 
canting  angles.  This  result  is  equivalent  to  that  corresponding  to  a  medium  having 
a  fraction  p  ^  of  seatterers  w ith  a  fixed  orientation  angle  or  and  a  fraction  1  -ft  (lf  with 
a  uniform  random  distribution  of  orientation  angles.  If  p  and  <v  are  independent  of 
drop  size,  thev  can  be  factored  out  of  the  right  side  of  Eq.  (12).  The  quantities 
0  and  can  then  be  defined  as  functions  of  Doppler  velocity  by 


-  \~$ 
i'  e 


f  olr.ol  N(r)  (i  i  (r,0l  c‘i^<  r’  ^  dr 
0 


f  n(r,  o)  N(r)  G  dr 
0 


(14) 
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BACKSCATTER  AMPLITUDE  RATIO  •>•<* 


Figure  1.  Baekscatter  amplitude  ratio  yeJ  for  raindrops, 
incorporating  the  Spiihaus  size-shape  relation  and  based  on 
the  scattering  coefficients  calculated  by  Oguchi  and  Hosoya. 
Calculations  are  shown  for  0.5  mm  increments  of  equivalent 
spherical  diameter  for  wavelengths  of  7.  5  cm  and  8.  6  mm 
and  elevation  angles  of  0°  and  40°.  Rayleigh  scattering 
approximation  holds  for  D  <  5  mm  at  7.  5  cm  wavelength  and 
for  D  <  0.  6  mm  at  8.  6  mm.  Magnitude  of  ratio  is  approxi¬ 
mately  proportional  to  cos^  <p.  Dependence  of  the  phase 
shift  6  on  elevation  angle  is  negligible  for  Rayleigh 
scattering  but  significant  for  drops  with  D/A  near  0.  5. 


If  the  weighting  function  G  is  omitted  from  the  integrals,  then  the  definition  applies 
to  the  entire  array  of  hydrometeors.  The  use  of  r  and  if  as  separate  parameters 
is  of  uncertain  value,  as  the  two  cannot  easily  be  modeled  separately,  except  in 
the  case  of  Rayleigh  scattering  where  5  -  0  or  tr.  One  can,  however,  approxi¬ 
mately  separate  the  size  and  elevation  angle  dependences  of  the  quantities  in  Eq. 

(12),  since  6  is  approximately  independent  of  elevation  angle  and  t\(r,  0)  -  r^(r ,  0)  cos4  ij> 
where  r.(r,0)  is  the  amplitude  ratio  for  a  drop  viewed  perpendicular  to  its  axis  of 
rotational  symmetry,  that  is,  horizontally. 


1 


In  summary,  the  backscatter  amplitude  ratio  for  a  collection  of  scatterers 
results  only  from  those  scatterers  having  a  common  orientation.  The  amplitude 
ratio  for  those  scatterers  can  be  computed  either  from  the  scattering  amplitude 
coefficients  for  the  specified  radar  elevation  angle  or  (approximately)  from  the 
scattering  amplitude  coefficients  for  zero  elevation  angle.  We  use  the  former 
procedure  in  our  computations 

2.3  Backscatter  Power  Ratio 

We  follow  the  approach  of  Warner  and  Rogers**  in  formulating  the  backscatter 
power  ratio.  This  quantity  is  the  sum  of  two  terms,  corresponding  to  the  preferen¬ 
tially  oriented  fraction  and  the  randomly  oriented  fraction  (1  -p  >,  respectively. 
The  power  ratio  can  be  expressed  as 

t-,2(r.(S)  -  ^(r)  cos 1  0  t-  [  1-p  ^  (r>|  (8M5iJ  ;'n>,0)  (15) 

where  0>  is  the  backscatter  power  r  atio  for  a  drop  viewed  at  horizontal  incidence 

1  •> 

and  the  relation  v .(r,  o)  •-  v.h-.O)  cos'  o  is  assur.  ed.  The  Ido  'or  8' 15  results  from 
4  1 

integrating  cos  ip  over  all  orientations.  Alternately,  we  can  calculate  the  quantity 
2 

rAr.o)  direct l v  from  the  scattering  amplitude  coefficients.  Using  this  procedure, 
we  compute  the  power  ratio  bv  the  formula 

t'"(r,ot  --  i>  (>(r)  u?  (r.  0)  *•  [l-p(y(r)]  (8/ 15)  r“’(r,  0)  (16) 

2 

The  power  ratio  r.(r.O)  calculated  from  the  scattering  amplitude  coefficients  of 
Oguchi  and  Hosoya  is  shown  in  Figure  2  for  two  wavelengths  as  a  function  of  fall 
speed  V(r)  which  was  calculated  from  the  formula  of  Best.  While  one  could  not 
obtain  this  curve  from  actual  observations  at  zero  elevation  angle,  because  the 
Doppler  velocity  component  due  to  fall  speed  is  proportional  to  sin  o,  one  could 
obtain  an  approximation  to  the  curve  by  plotting  (Sj/S^l/cos  o  against  v/sin0.  The 
result  would  be  identical  to  Figure  2  in  the  ideal  case  of  a  completely  oriented 
medium  with  zero  air  velocity.  Non-zero  air  velocity  would  result  in  a  displace¬ 
ment  along  the  abscissa.  Propagation  effects  and  air  velocity  variance  result  in 
changes  which  are  discussed  and  illustrated  in  Section  3. 


12.  Best,  A.  C.  (1950)  Empirical  formulae  for  the  terminal  velocity  of  water  drops 
falling  through  the  atmosphere,  Quart.  J.  Roy.  Meteor.  Soc.  76:302-311. 
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FALL  SPEED  (m  *c-’l 


Figure  2.  Scattering  power  ratio  at  zero 
elevation  angle  as  a  function  of  raindrop 
fall  speed  for  radar  wavelengths  of  7.5  cm 
and  8.  6  mm,  based  on  the  calculations  shown 
in  f  igure  1.  Power  ratio  for  8.8  mm  shows 
sharp  increase  near  8.  8  m  sec'^  fall  speed, 
corresponding  to  drops  of  diameter  near  half 
the  wavelength. 


The  foregoing  results  allow  the  spectral  functions  to  be  expressed  in  terms  of 


77, . 


the  subsidiary  functions  r  (v),  P<v),  and  S’fv),  so  that  Fqs.  (9)  and  (10)  become 

-  (^(v)  +  4p2  +  4p  ^  Re[po  j(Xi  2t)  F(v)  e|(^(v>i  2  (17) 

S../S.  -  n  T'(v)  el(5(v,i2«)  X  2neJ(X*2T>  . 


(18) 


2.4  Computational  Parameters  and  Procedure* 


Our  model  requires  as  input  parameters  the  rainfall  rate,  air  velocity  standard 
deviations,  radar  elevation  angle  0,  preferential  orientation  factor  pfl,  and 
propagation  distance.  We  use  the  scattering  amplitude  coefficients  and  rates  of 
differential  attenuation  (dB  km  *)  and  differential  phase  shift  (deg  km  *)  from 
Oguchi  and  Hosoya.  Consequently  the  model  computations  are  restricted  to  the 
frequencies  (4,  6,  11,  19.3,  30,  and  34.  8  GHz)  and  the  elevation  angles  (0*,  20®, 
10°,  and  60°)  for  which  the  amplitude  coefficients  were  calculated  and  the  rainfall 
rates  for  which  A.4  and  i$wero  calculated  (0.  25,  1.  25,  2.  5,  12.5,  25,  50,  100,  and 
150  mm  hr"*).  We  use  the  Marshall  and  Palmer  raindrop  size  distribution  func¬ 
tion  and  the  fall  speed  formula  of  Best.  For  simplicity  we  assume  at  =  7=0. 

We  calculate  the  term  pr1*  for  the  specified  propagation  distance.  We  then 
calculate  three  integral  quantities  at  intervals  of  0.  1  m  sec  *  between  -2  and 
!  6  m  sec  *; 


r 

max 

I.(v'  -  /  u(r,  o>  N(r)  G(Y(r)  sin  o,  li )  dr 

r  . 
min 


09) 


i2<v>  -  7 


u(r,0>  N(r)G(V(r)  sin  o. 


o  o 

L)  l'  (  r,  O)  +  (1-f)  )  -r?-  V 

t  a  i  y  a  15 


f(r.0») 


dr 


(20) 


max 


j  5.  (r,  0) 


/  .I  11  ,  VJI 

a(r,0)  N(r)  G(V(r)  sin  o,  -)  P.  (r,  q>)  e  dr  . 


(21) 


_  3 

We  use  r  =  0.  1  mm,  r  =3.25  mm,  and  Ar  =  3.  15  X  10  mm,  interpolating 
min  max  & 

the  scattering  coefficients  linearly  between  the  radii  at  which  they  were  calculated. 
The  power  spectra  and  cross -spectrum  are  then  calculated  as  follows; 


Si<y)  -  I  <v>  +  4p2I  <v)+4Re  [pe  *xI2(v)] 


(22) 


S2(v)  =  Ij(v)  (23) 

Sj2<v)  ■  Ig(v)  +  2pe^XIj(v)  .  (24) 
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From  these  we  form  the  secondary  spectral  functions  Sj/Sg  and  We  inte¬ 

grate  the  three  primary  spectral  functions  to  obtain  an  approximation  of  the  total 
reflectivity  as  defined  in  Eq.  (8).  the  circular  depolarization  ratio,  the  cross¬ 
covariance  and  cross -correlation,  and  the  Doppler  velocity  mean  and  standard 
deviation  in  each  channel. 

Noise  in  the  spectral  functions  can  be  simulated  in  two  ways.  Receiver  noise 
is  simulated  by  adding  to  both  power  spectra  a  uniform  noise  spectral  density  at  a 
specified  decibel  level  below  the  peak  of  the  spectrum  Sg-  This  corresponds  to 
radar  system  noise  with  a  Delta  function  autocorrelation  and  zero  cross -correlation 
between  the  noise  and  the  signal.  Because  this  noise  in  either  channel  is  assumed 
to  be  uncorrelated  with  the  signal  and  with  the  noise  in  the  other  channel,  it  has  no 
effect  on  the  cross -spectrum.  Noise  due  to  spectral  leakage,  as  would  be  encoun¬ 
tered  in  spectrum  analysis  procedures  applied  to  experimental  data,  is  simulated 
by  adding  to  each  function  a  uniform  noise  spectral  density  at  a  specified  decibel 
level  below  the  peak  magnitude  of  that  function.  In  the  case  of  the  cross -spectrum , 
the  constant  is  complex  and  has  the  same  argument  (or  phase)  as  the  peak  of  the 
cross  -sped  rum .  In  the  analysis  of  experimental  data  the  noise  spectral  density 
is  likely  to  he  dominated  more  often  by  spectral  leakage  than  by  receiver  noise. 

The  most  common  situation,  however,  is  likelv  to  combine  the  effects  of  spectral 
leakage  with  radar  svstem  noise  which  is  correlated  between  the  two  channels  due 
to  its  origin,  for  example,  in  a  local  oscillator. 

The  computations  illustrated  in  the  following  section  were  performed  at  Air 
Force  Geophysics  Laboratory  on  a  Perkin -Kil  ter  HU42  computer  and  the  results 
were  plotted  hv  a  California  Computer  Products  CALCOMP  1055  plotter.  Some 
apparently  anomalous  features  of  the  plotted  functions  are  due  to  the  numerical 
limits  of  the  computer,  that  is,  at  the  edges  of  the  Doppler  velocity  domain  of  the 
functions.  Numerical  accuracy  of  the  results  is  also  affected  bv  the  procedures 
we  used,  including  the  trapezoid -rule  integration,  the  radius  limits,  and  the*  linear 
interpolation  of  the  scattering  oeffieients. 

3.  RKSl'I.TS 

The  initial  objective  of  the  numerical  simulation  was  to  duplicate  the  theoretical 
forms  of  the  amplitude  and  power  ratios  (Figures  1  and  2),  determining  the  values 
of  the  computational  parameters  for  which  this  can  best  be  done  and  the  ways  in 
which  changes  in  computational  parameters  afreet  the  forms  of  the  output.  Compu¬ 
tations  were  performed  for  10’  elevation  angle,  2.5  and  25  mm  hr  rainfall  rates, 
and  several  values  of  air  velocity  variance,  propagation  distance,  and  noise.  This 
elevation  angle  was  chosen  to  yield  the  best  compromise  between  variance  due  to 


t 


2 

Doppler  fall  speed,  varying  as  sin  0,  and  backscatter  amplitude  ratio,  varying  as 
2 

cos  0.  Most  of  the  calculations  were  done  for  8.  6  mm  wavelength. 

Some  of  the  integrated  parameters  are  shown  in  Table  1.  The  circular  de¬ 
polarization  ratio  (CDR)  and  Doppler  mean  velocities  are  smaller  for  the  lower 
rainfall  rate,  due  to  the  greater  proportion  of  smaller,  more  spherical  drops. 
Differential  attenuation  and  differential  phase  shift  are  of  comparable  importance 
at  this  wavelength,  as  is  evident  in  the  vector  change  of  the  cross -covariance  ampli¬ 
tude  ratio  (CCAR)  with  propagation  distance.  Because  of  the  relative  phases  of  the 
scattering  and  propagation  terms,  the  initial  effect  of  propagation  is  to  decrease 
the  magnitude  of  the  spectrum  Sj,  especially  at  the  low -velocity  end.  Hence  at 
short  propagation  distances  the  CDR  decreases  and  the  velocity  increases.  At 
greater  propagation  distances,  the  CDR  increases  and  the  velocity  v  ^  asymptotically 
approaches  v  .  The  asymptotic  behavior  is  more  rapid  at  25  mm  hr’*  rainfall  rate 
because  of  the  larger  differential  attenuation  and  phase  shirt. 


Table  1.  Circular  Depolar'- zation  Ratio  (CDK),  Cross-Covariance  Amplitude 
Ratio  (CCAR),  and  lloppler  Mean  Velocities  as  1- unctions  of  Rainfall  Rate  (R> 
and  Propagation  Distance  (d),  for  8,(5  mm  Wavelength  and  40®  elevation  Angle 


R 

d 

■ 

GPR 

CCAR 

V1 

V2 

(mm  hr  S 

(km ) 

MR) 

magnitude 

n 

phase 

t,-.  sec"1) 

,  P 

(m  sec  ) 

(abs.  dute ) 

(dll) 

(dee ) 

2.  5 

0 

-28.  80 

0.  0342 

-14.  66 

147. 0 

4. 06 

3.70 

0.  5 

-28.  93 

0.  0337 

-14.77. 

140.  7 

4.  09 

1. 0 

-28.  97 

0.0335 

14.5  5 

134.  3 

4.  12 

2.  0 

-28.75 

0.0345 

-14.  62 

121.7 

4.  14 

4.0 

-27. 44 

0.0407 

-13.90 

100.  6 

4.  09 

25 

0 

-27.  50 

U. 0306 

- 14.  02 

137. 0 

4.  5  1 

4.  26 

0.  5 

-27. 55 

0. 0303 

-  T'..  05 

80.  5 

4.  64 

1.  0 

-23.  44 

0.  0657 

-11.33 

50.  3 

4.  46 

2.  0 

-17.38 

0.  135 

-8.70 

33.  3 

4.33 

_ 

The  following  subsections  discuss  the  effects  of  turbulence,  differential  propa¬ 
gation,  and  noise  in  more  detail.  Attention  is  placed  primarily  on  the  possibilities 
of  separating  Doppler  velocity  components  due  to  air  velocity  and  fall  speed,  sepa¬ 
rating  scattering  and  propagation  effects,  and  interpreting  the  spectral  functions  in 


the  presence  of  noise.  The  final  subsection  discusses  and  illustrates  these  effects 
in  relation  to  Rayleigh  scattering. 


3.1  Turbulence  Effect* 

Computations  for  scattering  at  8.  6  mm  wavelength  were  performed  with  several 
values  of  air  velocity  standard  deviation  to  illustrate  the  effects  of  air  velocity 
fluctuation  on  the  spectral  functions  and  to  attempt  to  generate  functions  approxi¬ 
mating  those  shown  in  Figures  1  and  2.  These  are  illustrated  in  Figures  3,  4. 
and  5  for  25  mm  hr  *  rainfall  rate  and  in  Figures  6  and  7  for  2.5  mm  hr  1  rainfall 
rate.  The  spectra  and  power  ratios  computed  with  2  =  0.  05  and  0.  1  show  varia¬ 
tions  which  are  probably  artifacts  of  the  linear  interpolation  of  the  scattering  coeffi¬ 
cients.  especially  between  5  and  6  m  sec  *  Doppler  velocity.  In  this  velocity 
region,  the  scattering  characteristics  change  rapidly  with  Doppler  velocity,  since 
the  raindrops  producing  these  effects  are  falling  within  7  0%  to  85%  of  the  maximum 
raindrop  fall  speed  and  the  Doppler  fall  speed  is  changing  slowly  with  drop  size. 

The  effects  of  the  numerical  procedure  with  small  values  of  2  are  more  strikingly 
illustrated  in  the  cross -spectral  amplitude  ratio  (CSAR)  shown  in  Figures  5a  and  7a. 
Because  of  the  rapid  variation  of  drop  size  with  fall  speed  for  drops  of  4  to  5  mm 
diameter,  the  full  rotation  of  the  amplitude  ratio  shown  in  Figure  1  could  only  be 
dupicated  if  the  spectral  functions  were  computed  at  much  smaller  velocity  inter- 
sals  with  correspondingly  smaller  values  of  2  and  with  non-linear  interpolation  of 
the  scattering  coefficients.  The  fact  that  this  rapid  rotation  of  the  amplitude  ratio 
vector  occurs  at  lower  doppler  velocities  than  would  be  expected  from  Figure  1  is 
difficult  to  explain;  this  is  probablv  related  to  the  weighting  of  the  cross -spectrum 
by  the  backscatter  cross-section  and  size  distribution  functions  and  to  the  irregu¬ 
larities  of  the  cross-section  functions  resulting  from  the  linear  interpolation.  At 
large  drop  sizes  the  amplitude  ratio  varies  less  rapidlv  with  drop  size,  and  the 
nature  of  the  computed  CSAR  near  6  m  sec  '  Doppler  velocity  in  Figures  5a  and  7a 
is  generally  preserved  with  higher  values  of  2  ,  although  it  appears  at  higher 
Doppler  velocities.  The  combined  effect  of  the  cross  -spectral  peak  (at  4.  5  m  sec  * 
Doppler  velocity  for  R  25  mm  hr  *1  and  the  portion  of  the  amplitude  ratio  curve 
for  the  largest  drops,  both  of  which  lie  in  the  upper  left  quadrant  of  the  complex 
plane,  is  to  shift  the  curve  of  the  CSAR  into  this  quadrant  with  increasing  2  and  leave 
only  a  relative  maximum  of  the  real  part  of  the  CSAR  as  a  remnant  of  the  rotation 
shown  in  Figure  1.  This  effect  is  evident,  for  example,  in  Figure  5c  for 
2  =  0. 5  m  sec  *  and  Figures  5e  and  7 e  for  2  -  1.0  m  sec  *.  As  2  increases  and 
the  cross -spectral  characteristics  are  spread  in  the  Doppler  velocity  domain,  the 
relative  maximum  appears  at  progressively  higher  Doppler  velocities.  Table  2 
documents  this  effect  and  illustrates  the  uncertainty  of  associating  the  Doppler 
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velocity  at  this  maximum  with  a  particular  drop  size,  that  is,  4.2  mm,  based  on 
Figure  1.  Such  an  association  of  a  particular  Doppler  velocity  with  a  Doppler  fall 
speed  component  would  be  a  key  to  extracting  the  Doppler  component  due  to  air 
Ve-i.eity. 


Table  2.  Doppler  Velocity  (m  sec  )  at  Maximum  of  Real 
Part  of  Cross -Speot  rat  Amplitude  Ratio,  for  8.  6  mm 
Wavelength  and  40’  Klevation  Angle 


R 

S  (m  sec  * ) 

(mm  hr  *) 

0.  05 

0.  10 

0.  5 

1.  0 

2,  5 

5.  6 

5.7 

7.2 

12.  3 

5.  6 

5.  6 

6.  2 

8.  7 

An  alternate  approach  to  'he  nn  blem  of  separating  Doppler  components  due  to 
fall  speed  and  air  velocity  is  to  find  a  Doppler  vt  Uk  itv  at  which  the  spectral  power 
ratio  (SPK>  or  the  CSA  R  does  not  change  with  11.  The  particular  value  of  the  SPR  is 
its  asymptotic  value,  namely,  the  CDR.  (A  comparable  procedure  i  an  not  be  applied 
•o  tlie  CSA  R  if  appreciable  uitferentiul  pi.  me  shift  in  scattering  is  present,  since  ttie 
asymptotic  value  of  the  CSAR,  that  is,  the  (.CAR,  does  not  generally  lie  on  the 
curve  defined  bv  the  CSAR.  )  The  Donpler  win,  ities  at  which  the  computed  SPR  is 
c'|ual  to  the  CDR  are  shown  in  Table  .3.  Con. par  mg  these  values  to  those  in  Table  1 

we  find  that  thev  are  within  l).  I  m  sec  1  of  the  Doppler  mean  velocity  v  for  each 

- 1  ^ 

rainfall  rate.  'I  hev  are  0.2  to  0.8  m  see  less  than  the  Doppler  velocities  of  the 
peaks  of  the  power  spectra  S,,.  If  one  takes  the  (  Dll  as  equal  to  the  power  ratio  of 
an  individual  drop,  then  the  CDRs  of  -28.80  and  -27.50  dB  irnplv  drops  of  1.8  and 
2.  1  mm  diameter,  with  Doppler  fall  speeds  of  3.  87  and  4.27  m  sec  ,  respectively. 
The  difference  of  these  values  and  the  respective  v.,  is  the  estimate  of  Doppler  air 
velocity,  which  is  zero  in  this  simulation.  Further  computations  will  determine  the 
a.  curacv  with  which  this  technique  van  he  applied  at  different  rainfall  rates.  Since 
both  the  CDR  and  v^  are  ratios  of  integrated  spectral  functions,  neither  is  affected 
bv  variations  in  N^,  the  coefficient  of  the  exponential  drop  size  distribution.  De¬ 
viations  from  this  assumed  distribution  may  affect  the  results,  as  will  variations 

13 

in  average  drop  shape  due  to  oscillation.  Recent  research  bv  Beard  et  al  indi¬ 
cates  that  the  mean  axial  ratios  of  oscillating  drops  arc  larger  than  the  equilibrium 

9 

values  derived  by  Pruppacher  and  Fitter,  which  in  turn  are  larger  than  the  values 


13.  Beard,  K.  V. ,  Johnson,  D.  B. ,  and  Jameson,  A.  R.  (1983)  Collisional  forcing 
of  raindrop  oscillations,  J,  Atmos,  Sci.  40:455-46  2. 
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derived  by  Spilhaus,  which  were  used  by  Oguchi  and  Hosoya  and  in  the  present 

model.  As  detailed  quantitative  results  of  the  drop  oscillation  studies  become 

available,  it  should  be  possible  to  revise  the  spectral  model  accordingly. 


Table  3.  Doppler  Velocity  (m  sec  )  at  Which  Spectral 
Power  Ratio  Equals  Circular  Depolarization  Ratio,  for 
8.  6  mm  Wavelength  and  40°  Elevation  Angle 


R 

2  (m  sec  * 

(mm  hr  * ) 

0.  05 

0.  10 

0.  5 

1.0 

2.  5 

3.  8 

3.  8 

3.7 

3.  6 

25 

4.  3 

4.  3 

4.  2 

4.  2 

3.2  Propagation  Effects 

The  effect  of  propagation  through  rain  is  to  change  the  relative  magnitude  of 
the  signal  in  the  transmission  channel  bv  a  component  which  has  a  range-dependent 
piiase  and  amplitude  relation  to  the  signal  in  the  opposite  channel.  This  effect  is 
most  distinct  in  the  cross -spectral  amplitude  ratio,  which  is  displaced  in  the  com¬ 
plex  plane  by  an  amount  equal  to  the  propagation  term  2peJ\  Examples  of  the 
spectral  power  ratio  and  the  cross -spectral  amplitude  ratio  in  the  presence  of 
propagation  effects  arc  shown  in  Figures  8,  9,  10,  and  11.  (These  figures  also 
illustrate  effects  of  noise,  to  be  discussed  below.  )  The  determination  of  the  propa¬ 
gation  term  in  the  received  signals  is  of  importance  in  that  it  permits  a  characteri¬ 
zation  of  the  propagation  medium,  either  as  a  whole  or  between  range  samples,  and 
permits  the  detailed  quantitative  interpretation  of  baekscatter  phenomena  in  signals 
that  are  "contaminated"  bv  propagation  effects. 

At  8.  6  mm  wavelength  the  SPR  is  not  usable  for  this  purpose,  since  one  does 

not  know  a  priori  the  relative  phase  angles  of  the  baekscatter  and  propagation  terms, 

that  is.  6  and  \.  (The  case  at  10  cm  wavelength  is  simpler,  since  5»  v  and 

\  -  tr/2.  )  Therefore  the  nropagation  term  must  be  deduced  from  the  CSAR.  In 

principle,  one  can  take  the  value  of  the  CSAR  at  the  lowest  Doppler  velocity  as  an 

estimate  of  the  propagation  term,  on  the  assumption  that  the  propagation  term 

dominates  the  baekscatter  term  in  that  region.  This  method  yields  an  imprecise 

result  in  the  presence  of  noise,  since  the  baekscatter  amplitude  ratio  does  not  reach 

zero  magnitude  within  the  domain  of  detection  of  the  power  spectrum  S2.  Thus  one 

must  establish  a  minimum  velocity  (v  .  )  for  the  detectability  of  the  spectrum 

nun  r  2 

at  which  one  takes  the  value  of  the  CSAR  as  the  estimate  of  the  propagation  term 
2pe^\  From  Figures  9  and  11  one  can  see  that  the  most  consistent  estimate  of  the 


propagation  term  can  be  obtained  from  the  value  of  the  CSAR  at  the  point  where  the 

real  part  is  a  maximum  near  the  lower  end  of  the  Doppler  velocity  domain.  Using 

a  fixed  radar  noise  level,  as  in  Figures  9a  and  11a,  one  finds  that  this  point  on  the 

curve  corresponds  to  a  value  of  S,  approximately  10  dB  above  the  noise  level.  At 
- 1  £ 

25  mm  hr  rainfall  rate  and  with  a  radar  noise  level  50  dB  below  the  peak  of  S„, 

.  1  i 
as  in  Figure  9a,  v  -  -  0.  2  m  sec  ;  with  a  radar  noise  level  40  dB  below  the 
6  min  _j 

peak  of  Sg,  v  .  -  0.  5  m  sec  .  Values  of  CSAR  at  the  minimum  velocity  for 

detectability  with  -40  dB  radar  noise  are  tabulated  in  Table  4,  where  they  can  be 
compared  to  the  corresponding  values  of  the  propagation  term.  The  estimates  of 
the  propagation  term  have  an  absolute  error  of  about  0.  007  5  independent  of  range 
for  R  2.  5  mm  hr  but  more  than  0.01  and  increasing  with  range  at 
R  -  25  mm  hr  \  The  error  for  the  case  R  -  2.  5  mm  hr  *  corresponds  to  one-way 
differential  atte-'t;  Rion  and  phase  shift  of  0.  065  dB  and  0.  038  deg.  While  these  may 
be  acceptably  small  for  some  purposes,  they  can  introduce  unacceptably  large  error 
in  the  resulting  "stimate  of  the  scattering  term  p  i'"(v)e^^v\  Comparison  of  cal¬ 
culations  of  the  CSAR  for  various  propagation  scenarios  shows  that  for  d  >  1.5  t.m 
at  2.  5  mm  hr  1  and  for  1  >  0.  2  km  at  25  mm  hr  *  the  "corrected"  CSAR  is  a  better 
approximation  to  the  propagation -free  CSAR  than  is  the  "uncorroctod"  CSAR.  In 
other  words,  the  correction  of  the  CSAR  for  propagation  effects  is  best  when  t he 
error  of  estimating  the  propagation  term  is  !--ss  than  the  propagation  coefficient  2p. 
Further  analysis  will  lit'  required  to  determine  the  implications  of  these  results 
for  the  derivation  oi  the  "corrected"  SPR  as  an  estimate  of  the  propagation-free 
SPK,  the  fraction  of  oriented  scatterers  .  and  other  barkscatter  parameters 
Use  of  the  method  described  with  empirical  data  should  yield  results  com¬ 
parable  to  those  shown  above  for  observations  of  a  backscattenng  medium  in  which 
the  smallest  scatterers  are  spherical.  In  general,  if  the  magnitude  oT  the  propa¬ 
gation-free  CSAR  were  zero  in  any  Doppler  velocity  band,  due  to  the  presence  of 
either  spherical  scatterers  or  randomly  oriented  scatterers,  one  could  use  the 
actual  value  of  the  CSAR  there  as  a  measure  of  the  propagation  term.  In  such  a 
case,  however,  one  would  ne  ed  an  independent  determination  of  such  a  feature  in 
the  propagation-free  CSAR,  perhaps  by  an  observation  at  very  close  range.  Such 
an  observation  might  be  obtained,  for  example,  in  snow  but  not  in  a  "bright  band" 
situation.  In  snow,  however,  one  could  not  be  sure  that  the  physical  parameters 
deduced  near  the  ground  would  be  representative  of  those  at  higher  altitude. 
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Table  4.  Lower  Doppler  Velocity  Limit  of  Detection  of  Power  Spectrum  S2  With 
-40  dB  Noise  Level  and  Corresponding  Value  of  Cross -Spectral  Amplitude  Ratio 
(CSAR),  for  8,  6  mm  Wavelength  and  40°  Elevation  Angle 


R 

(mm  hr  *) 

CSAR(v  .  ) 
min 

2pej* 

2.  5 

0 

-0. 2 

-0.007  51  +  j  0.00067 

0 

0.  5 

-0.  00513  +  j  0.  00307 

0.00265  4  j  0.00267 

1.  0 

-0.  0027  6  +  j  0.00546 

0.  00530  +  j  0.  00534 

2.  0 

0.  00199  +  j  0.  0103 

0.  0106  4  j  0.  0107 

25 

0 

a  5 

-0.  0106  +  j  0.  00122 

0 

mm 

0.  0219  +  j  0.  0120 

0.0355  +  j  0.0118 

D 

0.0545  *  j  0.0228 

0.  0709  4  j  0.0235 

8.3  Noise 

Effects  of  radar  system  noise  on  the  spectral  power  ratio  are  shown  in  Figures 
4  and  6.  In  the  presence  of  such  noise,  the  SPR  is  an  accurate  representation  of 
its  noise-free  form  only  in  the  velocity  region  in  which  the  power  spectrum  in  the 
transmission  channel,  Sj,  is  more  than  about  10  dB  above  the  noise  level.  Thus, 
for  the  cases  illustrated  here,  the  SPR  is  of  little  value  if  the  radar  noise  level 
were  higher  than  about  50  dB  below  the  peak  of  S,;.  Because  the  radar  system  noise 
is  assumed  to  be  uncorrelated  between  the  two  polarization  channels,  it  does  not 
affect  the  cross -spectrum ,  and  the  CSAR  is  therefore  definable  across  the  entire 
Doppler  domain  within  which  S2  is  above  the  noise  level.  Radar  noise  effects  can 
be  seen  by  comparing  f  igures  5c  and  5d,  Figures  5e  and  5f,  Figures  7c  and  7d, 
and  Figures  7o  and  7f.  It  appears  that  in  the  presence  of  radar  noise  the  CSAR 
is  u  more  useful  function  for  analysis  than  is  the  SIR.  For  example, 
although  the  relative  maximum  of  the  real  part  is  lost  when  -  1.0  and  the  radar 
noise  level  is  at  -50  dB  (for  example.  Figure  5f).  the  form  of  the  function  is  pre¬ 
served  aceuratelv  for  Doppler  velocities  up  to  about  8  m  sec  .  The  quantitative 
details  of  the  CSAR  are,  of  course,  dependent  on  the  orientation  parameter  p . 

If  one  simulates  the  noise  as  being  due  to  spectral  leakage,  as  described  in 

Section  2,  then  the  SPR  and  CSAR  approach  values  of  S,  / S„  and 

1  max  2  max 

S,n  /S„  at  the  extremes  of  the  Doppler  velocity  domain.  These  effects 
12  max  2  max 

are  illustrated  in  Figures  8b  and  10b  for  the  SPR  and  Figures  9b  and  lib  for  the 
CSAR.  This  simulation  permits  the  definition  of  Sj  and,  hence,  the  SPR  over  a 
wider  velocity  domain  than  in  the  ease  of  simulated  radar  noise  described  above. 
Because  the  spectral  noise  affects  the  cross -spectrum,  the  CSAR  is  definable  over 
a  narrower  velocity  domain  than  in  ihe  case  of  simulated  radar  noise. 
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In  practice,  one  will  be  interested  in  situations  for  which  the  CDR  is  as  low  as 

.3 

10  ,  that  is,  -30  dB,  assuming  that  the  radar  antenna  permits  such  measurements. 

In  order  to  have  useful  signal  in  the  transmission  channel,  Sj  must  be  more  than 
10  dB  above  the  radar  noise  level  within  an  appreciable  velocity  domain.  The 
spectrum  in  the  opposite  channel,  having  a  peak  value  of  about  30  dB  higher  than  the 
peak  of  Sj.  will  be  more  than  40  dB  above  the  radar  noise  level  within  a  similar 
velocity  domain.  Since  common  data  windows  used  in  conjunction  with  fast  Fourier 
transforms  for  spectral  estimation  will  suppress  spectral  leakage  to  about  40  or 
45  dB  in  frequency  side  lobes,  the  domain  of  the  power  spectrum  will  generally 
be  limited  by  spectral  noise  rather  than  bv  radar  system  noise.  If  the  spectrum  Sj 
were  also  limited  bv  spectral  noise  at  the  40  dB  l<  vel,  then  the  radar  noise  level 
would  have  to  be  more  than  70  dB  below  the  peal;  of  S.,.  implying  the  possibility  of 
receiver  saturation  in  the  channel  opposite  to  the  transmission  channel.  Hence  Sj 
is  most  like'v  to  be  limited  by  radar  system  noise.  If  the  radar  is  lesigned  with 
high  isolation  between  chan  ids,  then  Hie  cross -spo  t rum  will  alwavs  be  limited  bv 
spectral  leakage  noise. 


3.4  Rawli-igli  Scatterin'! 


Calculations  analogous  :■>  the  foregoing,  mil  for  a  wave-length  .  5  in,  art- 
illustrated  in  Figures  1C- 1 H .  At  this  wavelength  the  scattering  di IT-  renltal  phase 
shift  6  is  near  r,  for  equivalent  drop  li'imet.-rs  of  4  nun  anti  loss.  These  sizes 
generally  dominate  the  computed  .spectral  functions,  especially  for  moderate  values 
of  il,  so  that  th.*  scattering  term  in  the  cross -spec1  ral  amplitude  ratio  is  pre¬ 
dominantly  real.  Attenuation  anti  differential  attenuation  are  small,  so  that  the 

j  V 

propagation  term  2pi  is  due  almost  entirely  to  differential  phase  shift.  The  angle 
\  ranges  from  0.48  if(86.8*)  to  0.  4!>~  (88.  9  ‘  1  for  rainfall  rates  of  150  to  0.25  nun  hr 
At  25  mm  hr  ^  rainfall  rate  and  40°  eievatimi  angle  with  -  w  0.5  m  set  *  (Figure  14b) 
the  angle  Z  lies  between  179.  5*  anti  -174..*'  a<  toss  the  velocity  Interval  0.3  to 
7.7  m.  see  ,  where  the  spectrum  is  mere  than  5  dB  above  a  -40  dB  radar  noise 


level. 

If  the  real  part  of  the  CSAK  is  interpreted  as  due  to  scattering  and  the  imaginary 
nart  due  to  propagation,  the  error  of  estimating  the  propagation  term  will  be  given 
bv  the  quantity  [Re"(2pe-1')  •  Im  (p  (v  .  I  e1 '  min  )]  .  At  a  propagation 

distance  of  2  km  through  heavy  (25  mm  hr  )  rain,  illustrated  in  Figure  16,  the 
magnitude  of  error  is  about  0.0016,  due  largely  to  Hie  imaginary  part  of  the 
scattering  term.  The  simpler  relationship  of  the  scattering  and  propagation  terms 
at  the  longer  wavelength  should  facilitate  the  identification  of  non-Rayleigh  scatterers 
and  the  presence  of  non-zero  canting  angles  in  the  scattering  and  propagation  media. 
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1  igure  5a.  Cross-spectral  amplitude  ratio  for  40° 

elevation  angle  in  25  mm  hr  1  rain  with  0.  05  m  sec 
and  no  noise.  Segmented  appearance  of  function  beyond 

about  5.  5  m  sec  *  Doppler  velocity  is  due  to  rapid 
variation  of  drop  size  and  shape  with  Doppler  fall  speed 
in  this  region. 


Figure  5c.  Same  as  Figure  5a 

but  with  £  =  0.  5  m  sec  1 
Smoothing  due  to  increasing 
width  of  velocity  weighting 
function  is  evident. 


Figure  5d.  Same  as  Figure  5c 
but  with  radar  noise  level  of 
50  dB  below  the  peak  value 


Figure  7c.  Same  as 
Figure  7a  but  with 

S  =  0.  5  m  sec  * . 


0.02 


Figure  7d.  Same  as  Figure  7c 
but  with  radar  noise  level  of 
50  dB  below  the  peak  value 
of  Sj,.  Disappearance  of  the 

relative  maximum  of  the  real 

part  near  7  m  sec  1  is  due  to 
the  narrower  Doppler  velocity 
domain  of  the  spectral  functions 
at  this  rainfall  rate  compared 

to  those  at  25  mm  hr  .  Compare 
this  figure  to  Figure  5d. 


Figure  7e.  Same  as 
Figure  7a  but  with 

=  1.0  m  sec 


Figure  7f.  Same  as  Figure  7e 
but  with  noise  level  of  50  dB 
below  the  peak  value  of  S? 


Figure  8a.  Spectral  power 
ratio  for  40°  elevation  angle 

in  25  mm  hr  rain  with 

^  =  0. 5  m  sec  ,  radar  noise 
level  of  50  dB  below  the  peak 
value  of  $2.  and  propagation 

distances  of  0,  0.  5,  and  1  km. 
Changes  with  distance  are  due 
to  the  relative  magnitudes  and 
phases  of  the  backscattering 
and  propagation  terms. 


Figure  8b.  Same  as  Figure  8a 
with  no  radar  noise  but  with 
spectral  noise  in  each  spectrum 
at  40  dll  below  the  respective 
peak  value.  Useful  domain  of 
the  function  is  increased  relative 
to  that  in  the  presence  of  radar 
noise. 


Re(Sw/S,' 


Figure  9a.  Cross-spectral 
amplitude  ratio  for  40° 
elevation  angle  in 

25  mm  hr  1  rain  with 

=  0.  5  m  sec  radar 
noise  level  of  50  dB 
below  the  peak  value  of  Sg 

and  propagation  distances 
of  0,  0.  5,  and  1  km.  Since 
p  =  1,  the  squared 

magnitude  of  this  function 
is  approximately  equal  to 
the  spectral  power  ratio 
shown  in  Figure  8a  for  the 
same  values  of  propagation 
distance. 


Figure  9b.  Same  as 
Figure  9a  with  no  radar 
noise  but  with  spectral 
noise  at  40  dB  below  the 
peak  values  of  and 

of  S^.  Compare  this 
figure  with  Figure  8b. 
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Figure  ]0a.  Spectral  power 
ratio  for  40“  elevation  angle 

in  2.  5  mm  hr  *  rain  with 

w  =  0.  5  m  sec  ,  radar  noise 
level  of  50  dB  below  the  peak 
value  of  Sg.  and  propagation 

distances  of  0,  2,  and  4  km. 
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Figure  10b.  Same  as  Figure  10a 
with  no  radar  noise  but  with 
spectral  noise  in  each  spectrum 
at  40  dB  below  the  respective 
peak  value. 
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Figure  11a.  Cross-spectral 
amplitude  ratio  for  40°  elevation 

angle  in  2.  5  mm  hr  1  rain  with 

2  =  0.  5  m  sec  1 ,  radar  noise  level 
of  50  dB  below  the  peak  value  of 
Sg.  and  propagation  distances  of 

0,  2,  and  4  km.  Compare  this  figure 
to  Figure  10a  to  see  the  roles  of  the 
relative  magnitude  and  phase  of  the 
backscattering  and  propagation  terms 
in  defining  the  range-wise  changes 
of  the  spectral  power  ratio. 


Figure  lib.  Same  as  Figure  11a 
with  no  radar  noise  but  with 
spectral  noise  at  40  dB  below  the 
peak  values  of  S^  and 

Compare  this  figure  to  Figure  10b. 
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Figure  12a.  Power  spectra  of 
signals  in  the  transmission  (1) 
and  orthogonal  (2)  channels  of 
a  circularly  polarized  7.5  cm 
wavelength  radar  at  40°  elevation 

angle  in  25  mm  hr  *  rain  with 

LI  =  0.  05  m  sec  .  Compare 
this  figure  to  Figure  3a,  which 
illustrates  non- Rayleigh 
scattering  with  identical  values 
of  all  parameters  except  the 
wavelength. 
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f  igure  12b.  Same  as  Figure  12a 
but  with  LI  =  0.  5  m  sec”  l. 
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Figure  14a.  Cross-spectral  amplitude  ratio  for  7.5  cm 
wavelength  and  40°  elevation  angle  in  25  mm  hr  1  rain  with 
v  =  0.  05  m  sec  '  and  no  noise.  Divergence  of  function 

from  negative  real  axis  beyond  about  5.  7  m  sec  '  Doppler 
velocity  is  due  to  non-  Kavleigh  scattering  effects. 


figure  14b.  Same  as  Figure  14a  but  with  11  =  0.  5  m  sec  . 


o 


Figure  14c.  Same  as  Figure  14b  but  with  radar  noise 
level  of  50  dB  below  the  peak  value  of  S£.  Effect  of 

noise  is  to  obscure  most  of  the  non-Kayleigh  scattering 

component  above  7  m  sec  *  Doppler  velocity. 
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Figure  15a.  Spectral  power 
ratio  for  7.  5  cm  wavelength 
and  40°  elevation  angle  in 

25  mm  hr  1  rain  with 

2  =  0.  5  m  sec  ' .  radar  noise 
level  of  50  dB  below  the  peak 
value  of  S2>  and  propagation 

distances  of  0,  2,  and  4  km. 
The  function  increases  mono- 
tonically  with  distance  because 
the  backscattering  and 
propagation  terms  differ  in 
phase  by  approximately  90°. 
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Figure  15b.  Same  as  Figure  15a 
with  no  radar  noise  but  with 
spectral  noise  in  each  spectrum 
at  40  dB  below  the  respective 
peak  value. 
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Figure  16a.  Cross-spectral  amplitude  ratio  for  7.  5  cm 
wavelength  and  40°  elevation  angle  in  25  mm  hr  1  rain 

with  2  =  0.  5  m  sec  \  radar  noise  level  of  50  dB  below 
the  peak  value  of  $2.  and  propagation  distances  of  0,  2, 

and  4  km.  Imaginary  part  of  this  function  taken  between 

1  and  5  m  sec  1  Uoppler  velocity  yields  an  excellent 
approximation  of  the  propagation  term,  which  is  due 
mainly  to  differential  phase  shift. 
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Figure  16b.  Same  as  Figure  16a  with  no  radar 
noise  but  with  spectral  noise  at  40  dB  below  the 
peak  values  of  S^  and  Sg. 


4.  CONCLUSIONS 


We  have  developed  a  model  for  simulating  spectral  functions  derivable  from  a 
coherent  polarization  diversity  radar  transmitting  circular  polarization  and  ob¬ 
serving  backscatter  from  rain.  The  spectral  functions  have  been  computed  for  radar 
wavelengths  of  8.  6  mm  and  7.  5  cm  and  for  several  values  of  air  velocity  variance 
and  propagation  distance.  We  performed  these  computations  for  40°  elevation 

angle,  for  jointly  maximizing  the  signal  amplitude  in  the  transmission  channel 
2 

(proportional  to  cos  0)  and  the  Doppler  velocity  variance  due  to  fall  speeds  (pro- 
2 

portional  to  sin  q>). 

The  results  illustrate  several  requirements  that  must  be  met  in  the  analysis 
of  this  type  of  radar  data.  The  spectrum  in  the  transmission  channel  must  be  more 
than  10  dB  above  the  noise  level  over  an  appreciable  velocity  domain  in  order  that 
the  spectral  power  ratio  be  useful  for  analysis.  The  cross-spectral  amplitude 
ratio  appears  to  have  more  analytical  value  because  it  is  not  as  severely  affected 

by  noise.  Large  values  of  air  velocity  variance,  for  example,  approaching 

2  -  2 

1  m  sec  ,  limit  the  usefulness  of  the  spectral  functions,  as  certain  features  that 
might  be  associated  with  particular  drop  sizes  are  thereby  shifted  in  Doppler 
velocity. 

A  fairly  accurate  estimate  of  the  Doppler  component  of  air  velocity  can  be  de¬ 
rived  at  8.  6  mm  wavelength  by  identifying  the  Doppler  fall  speed  component  at  v 2 
with  drops  of  size  corresponding  to  the  power  ratio  SjlVgl/SgfVg).  This  estimate 
is  dependent  on  the  size-shape  relation  assumed  in  the  model,  and  refinements  of 
the  specification  of  average  drop  shape  will  lead  to  model  results  more  closely 
applicable  to  experimental  measurements.  Use  of  this  procedure  in  the  presence 
of  differential  propagation  effects  will  lead  to  inaccuracies  unless  the  propagation 
term  is  estimated  accurately.  The  present  model  can  be  extended  to  permit  evalua¬ 
tion  of  the  resulting  uncertainties.  Calculations  at  7.5  cm,  although  not  as  com¬ 
prehensive  as  at  8.  6  mm,  indicate  that  comparable  results  can  be  obtained. 

Estimation  of  the  total  propagation  term  in  the  CSAK  yields  an  error  of  about 
0.  0075  to  0.  010  in  magnitude  at  8.  6  mm  wavelength,  which  may  not  be  acceptable 
for  some  purposes.  The  incremental  propagation  term,  that  is,  from  one  range 
gate  to  the  next,  can  be  evaluated  with  more  accuracy,  probably  exceeding  the 
accuracy  obtainable  from  the  CCAK  based  on  non-coherent  radar  data.  At  7.  5  cm 
wavelength,  where  the  scattering  amplitude  ratio  is  predominantly  real,  the  total 
propagation  term  can  be  estimated  with  much  more  accuracy  than  at  8.6  mm. 

Several  future  developments  and  applications  of  the  model  are  possible.  Fur¬ 
ther  calculations  can  be  made  with  the  existing  formulation  to  evaluate  uncertainties 
of  estimating  the  Doppler  air  velocity  component  and  the  propagation  term.  The 
noise  simulation  options  should  be  combined  to  more  closely  match  experimental 
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conditions.  A  signal  component  corresponding  to  radar  system  cross-polarization 
could  be  included.  Of  greatest  scientific  interest  (and  of  greatest  difficulty)  would 
be  the  modification  of  the  model  to  incorporate  the  newest  data  on  raindrop  shapes 
and  to  permit  the  simulation  of  mixed-phase  media. 
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